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Studies of the nature of cosmic ray particles at the highest energies are based on the measurement 
of extensive air showers. Most cosmic ray properties can therefore only be obtained from the 
interpretation of air shower data and are thus depending on predictions of hadronic interaction 
models at ultra-high energies. We discuss different scenarios of model extrapolations from accelerator 
data to air shower energies and investigate their impact on the corresponding air shower predictions. 
To explore the effect of different extrapolations by hadronic interaction models we developed an ad 
hoc model. This ad hoc model is based on the modification of the output of standard hadronic 
interaction event generators within the air shower simulation process and allows us to study the 
impact of changing interaction features on the air shower development. In a systematic study we 
demonstrate the resulting changes of important air shower observables and also discuss them in 
terms of the predictions of the Heitler model of air shower cascades. It is found that the results of 
our ad hoc modifications are, to a large extend, independent of the choice of the underlying hadronic 
interaction model. 



I. INTRODUCTION 

The nature of the highest energy cosmic ray particles 
is still an open question. Even with recent high expo- 
sure experiments delivering large quantities of data at 
ultra- high energies [lJ-Q, the fundamental problem re- 
mains to interpret the collected air shower data. This 
is mainly due to difficulties and limitations in modelling 
of hadronic interactions in air shower cascades. Whereas 
we can predict reliably electromagnetic interactions with 
QED, it is still not possible to calculate hadronic multi- 
particle production from first principles. Moreover, most 
of the interaction energies and phase space regions of rel- 
evance to the development of air showers are not directly 
accessible in accelerator experiments. Thus, phenomeno- 
logical models of hadronic interactions have to be used 
to extrapolate the available accelerator measurements to 
these unexplored phase space regions. 

Due to the lack of a theoretical framework that allows 
one to make quantitative predictions, the systematic un- 
certainties of these extrapolations are very difficult to 
estimate. This is a well known problem for, e.g. the 
determination of the systematic uncertainty of the pri- 
mary mass composition derived from air shower measure- 
ments Similarly, the reconstruction of the energy 
of the primary particle from detector data calibrated by 
Monte Carlo simulations is subject to largely unknown 
model uncertainties 

For estimating systematic uncertainties of air shower 
predictions, one first needs to know the uncertainties of 
the model predictions for multiparticle production, and 
secondly how these uncertainties propagate to predic- 
tions for different air shower observables. In this work 
we will address the latter by studying the relation of dif- 
ferent extrapolations of hadronic particle production to 
air shower predictions. 

We will briefly discuss the uncertainties associated 



with the modelling of hadronic interactions and argue 
that the differences of existing models do not cover the 
full range of the expected uncertainty at energies beyond 
the reach of colliders. To nevertheless explore a wide 
range of different extrapolations and to calculate the cor- 
responding air shower observables, we will introduce an 
ad hoc modification of existing models. The shower ob- 
servables under consideration are the depth of shower 
maximum and the number of muons at ground level-as 
the model-dependence of the interpretation of air shower 
data is mainly related to these two parameters (l5j. All 
other characteristics of air showers are closely related to 
those two features due to universality of the electromag- 
netic shower component (see [l6l - [2lT ] and Refs. therein). 
The results will also be compared to the expectations 
from the Heitler-Matthews cascade model 1231. 



II. UNCERTAINTIES IN MODELLING 
HADRONIC INTERACTIONS 

The development of extensive air showers (EAS) is very 
sensitive to the characteristics of hadronic interactions at 
ultra-high energies. Different hadronic interaction event 
generators are available for the simulation of EAS. 

Currently it is unclear whether the existing differences 
in model predictions for hadronic multiparticle produc- 
tion can be used to estimate the level of theoretical uncer- 
tainties related to our limited understanding of hadronic 
interactions in air showers. It is possible that existing 
model-differences are 

(A) larger than the actual systematic uncertainties. For 
example, progress in understanding hadronic par- 
ticle production and new data from accelerators 
allow us to update interaction models to obtain 
a more realistic description of particle production. 
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FIG. 1. Compilation of accelerator data of af£ t and B e \ [23]. The central line denotes the conventional extrapolation of these 
data to high energy. The upper and lower lines indicate a set of possible extreme extrapolations. In the left plot the conventional 
model is the soft pomeron parametrization by Donnachie and Landshoff [24|, while the lower curve is by Pancheri et al. [25| 
and the upper one is the two-pomeron model of Landshoff [261 |27| . The different scenarios in the right plot are from [2|| . 



Not all models arc updated regularly and the qual- 
ity of data description differs between the models. 

(B) smaller than the actual systematic uncertainties. 
The existing models are not covering the full phase 
space of possible interaction scenarios and param- 
eters. Moreover, new physics processes at higher 
energies, which are unknown now and thus miss- 
ing in current modelling approaches, could change 
extrapolations drastically. 

Frequently used models for the high-energy range are 
QGSJet II [pJML EPOS [HJ, and the somewhat older 
QGSJet 01 [Hill and Sibyll 2.1 Q. These mod- 
els are available in the air shower simulation package 
CONEX [HI that will be used for calculating the shower 
observables. Other models for hadronic interactions in- 
clude NEXUS [MHz!, HDPM [H|, DPMJET [H], and 
VENUS [4(J. These models are older or more limited in 
the scope of application and not considered here. 

Despite the different level of sophistication, the pre- 
dictions by Sibyll 2.1 and QGSJet 01 are not ob- 
jectively worse than those by QGSJetII and Epos, as 
many model aspects are assumptions that cannot be jus- 
tified by underlying fundamental theoretical constraints. 
Over the years model predictions and extrapolations have 
become more alike even though there is no theory for cal- 
culating e.g. cross sections from first principles [4l| . One 
has to be careful and should not consider this increas- 
ing similarity of model predictions as real convergence 
and significant decrease of the uncertainties. None of 
the models is able to consistently describe cosmic ray 
data (e.g. 0, H, IH, El])- In the energy range up to 
about 10 15 eV where various measurements on multipar- 



ticle production are available good tuning to many dif- 
ferent data sets should indeed lead to a convergence of 
the model predictions. However, at energies beyond that 
of collider experiments, the extrapolations can only be 
guided by theoretical end phenomcnological assumptions. 

In Fig. Q] accelerator data on the total proton-proton 
cross section er{^ and the elastic slope parameter B e \, 
defined by d<7 c i/dt| t=0 oc exp(— |t|j3 c i), are shown to- 
gether with different models that extrapolate these data 
to ultra-high energies. Converting these model extrapo- 
lations within the Glauber framework [HI, J4f| to proton- 
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FIG. 2. Uncertainty of the extrapolation of the proton-air 
cross section for particle production due to different models of 
the proton- pro ton cross section as calculated with the Glauber 
framework [28j ]. 
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air cross sections for particle production as needed for air 
shower simulation leads to a wide range of predictions at 
ultra- high energy [28J . This range is shown in Fig. as 
the shaded area. The extrapolations of aforementioned 
interaction models are also shown and it is obvious that 
they do not exhaust the possible range of uncertainties. 
In addition it can be seen that a measurement of the 
total cross section at LHC energy with an uncertainty 
smaller than a few percent has the potential to signifi- 
cantly reduce the uncertainties of the model extrapola- 
tions at cosmic ray energies. Furthermore, recent activ- 
ities of model development of the EPOS event generator 
clearly demonstrated that there exists a large freedom to 
change predictions at ultra- high energies 0, [47j ■ 

All this is demonstrating that the current situation is 
most likely such that one would underestimate the sys- 
tematic uncertainties of model extrapolations at ener- 
gies beyond the reach of accelerator experiments if one 
would just consider the extrapolation of existing event 
generators. This is supported by the fact that existing 
interaction models are not able to consistently describe 
cosmic ray observations. Therefore it is also not sur- 
prising that first LHC data [48l - |55j indicate some model 
deficiencies (56j . 

In the following we will consider different extrapola- 
tions of multiparticlc production to energies higher than 
10 15 eV, corresponding to a c.m.s. energy of about 2 TeV 
for proton-proton collisions. We will concentrate on gen- 
eral features of hadronic particle production that are 
most directly linked to air shower predictions. 



III. PREDICTIONS OF THE CASCADE MODEL 
FOR AIR SHOWERS 



n=1 




FIG. 3. Electromagnetic Heitler model, 
and their energy is 

E n = E /N n . (2) 

Defining the critical energy E° m - as the energy below 
which energy loss processes dominate over particle pro- 
duction (~ 85MeV in air), one can make the assumption 
that the shower maximum is reached when the energy 
of secondary particles reaches E° m \ Thus, two main 
shower observables are given by 

AW - E /E c c - m - ~ E and 

X max = X c n c ~ A c ln(E ) , (3) 

with n c = (ln2) _1 In (Eq/E^™-). Even this very simplis- 
tic model reproduces two important features of air show- 
ers: the number of particles at shower maximum A max 
is proportional to Eq and the depth of shower maximum 
A max depends logarithmically on the primary energy E . 



Simple cascade models, often referred to as Heitler 
models (57j . are providing some insight into how air 
shower observables are related to interaction physics on 
a microscopic level [22|, [5^, [59[ . The descriptive strength 
of Heitler modcls-dcspite their extreme simplicity-is re- 
markable. We will discuss the ability of such models to 
link the physics of interactions to air shower observables, 
but also point out the limitations of these models. 



A. Electromagnetic Heitler Model 

In the electromagnetic Heitler approximation only one 
particle type is considered and substitutes 7, e + and e~. 
It is assumed that an e.m. particle with energy E in- 
teracts after one splitting length A c = ln2Ao, where 
Xo ~ 37g/cm 2 is the e.m. radiation length, producing 
two secondaries with energy E/2, see Fig. [3] The num- 
ber of particles after each splitting length increases by a 
factor of two, and thus the number of particles at gener- 
ation n is 

A„ = 2" (1) 



B. Hadronic Extension of the Heitler Model 

The Heitler model can be extended to hadronic par- 
ticles by considering a cascade of pions interacting in 
air [22| . It is assumed that a charged pion (n + or 
7r~) of energy E interacts after one interaction length 
Ai (» 120g/cm 2 for pions of 10 - 1000 GcV and 
produces r n mu it charged pions and cn mu \t neutral pions, 
see Fig. [4] Here c = 1 — r is the pion charge-ratio and 
is in the Heitler framework typically defined to be 1/3. 
Neutral pions decay instantaneously into 27 and are thus 
lost from the hadronic cascade, while transferring energy 
into the electromagnetic cascade. The number of charged 
pions in the hadronic cascade at generation n is 

Af = (rn mult )" (4) 

and the energy per pion 

K = r-^V ■ (5) 

It is assumed that charged pions decay into muons as 
soon as their energy drops below E& ec . This is the 
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FIG. 4. Hadronic Heitler-Matthews model: pion cascade in 
air. 



characteristic energy where the interaction length be- 
comes larger than the decay length of pions in air, and 
is thus also depending on the atmospheric density pro- 
file. For example high up in the atmosphere at a depth 
of ~ 100g/cm 2 this energy is around -Edcc ~ 55 GcV 
and decreases to E^ ec ~ 18 GeV at sea level. In realistic 
air shower scenarios the energy E™ = Edec is typi cally 
reached within the range of 20 - 30 GcV [H [gO, EHf , 
which is then adopted to be the critical energy of pions 
E%, in analogy to the critical energy of electrons E^ m -. 
The number of generations required to reach E* = E^ is 



HEq/e?) 

lnn m „it. 



(6) 



and is for air showers in the range of 4 to 7 [62| . This 
yields a muon number from decaying charged pions of 



N ch = 

n c 
In rn„ 



lnr 



In n T 



In n r 



Eo 
E* 

■1, 



with 



(7) 



which depends on the secondary particle multiplicity and 
the pion charge ratio of hadronic interactions. Moreover, 
the electron number at shower maximum can be esti- 
mated from 



E e 



Eo 



Eft 



N„ 



(8) 



and is thus depending on the hadronic secondary multi- 
plicity and pion charge ratio via the muon number. To 
relate the electron number at the shower maximum to 
the electron number at observation level, N e , it is nec- 
essary to take the strong attenuation of the electromag- 
netic shower component after the shower maximum into 
account. In the limit where X D b s 3> l m ax the air shower 
size is in very good approximation exponentially atten- 
uated with the scale length A ~ 65g/cm 2 [glj]. We can 
thus use 



JV R ~JV B 



0/A 



which leads to 
ln7V e - 



— (A"obs — -^max)/A , 



(9) 



It is important to note that In N c depends only logarith- 
mically on JV e max but linearly from X max . It is shown in 
Section|V]that indeed the impact of the depth of X max on 
the electron number at ground level is dominating over 
a change following from Eq. ©. It is thus the strong 
longitudinal shower evolution that is responsible for the 
inability of the Heitler model to directly infer any depen- 
dence of N e on hadronic interaction characteristics. 

Within the hadronic Heitler framework it is not pos- 
sible to follow the parallel development of the hadronic 
and electromagnetic cascades. To estimate the position 
of the shower maximum one is forced to consider the elec- 
tromagnetic contribution from the first hadronic interac- 
tion only, which consists of 2 c n mu i t e.m. subshowcrs each 
of the energy -Eo/^muit- The shower maximum is then 



Ai + A e In 



En 



"mult E° 



(10) 



which is proportional to the hadronic interaction length 
and depends logarithmically on the multiplicity, but not 
on the charge ratio. This expression does not include 
the contribution to the electromagnetic cascade from the 
subsequent hadronic cascading process. However, the in- 
clusion of higher hadronic generations does not change 
the structure of Eq. (P~0|) . only the coefficients change 
(e.g. 0). 



C. Inelasticity in the Heitler Model 

In the Heitler model only equal energy particles of one 
type arc considered. This also excludes any account of 
leading particle effects or other secondary particle dis- 
tribution effects. To some limited extend it was accom- 
plished to incorporate the inelasticity of interactions in a 
Heitler-type cascade by Matthews |22|]. The inelasticity 
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(11) 



for X obs > X, 



is the fraction of the primary energy that is not car- 
ried away by the most energetic secondary particle, often 
referred to as leading particle. This energy fraction is 
available for the production of new secondary particles, 
mainly pions and kaons. The elasticity K e i = 1 — ftmci is 
the fraction of the primary energy that is carried by the 
leading particle. 

The main difference compared to the standard Heitler 
model is that after each interaction secondaries with two 
energy levels are generated: n mu it particles with energy 
of Kind Eq / n mu \t (of which c n mu it are neutral pions) and 
one particle with energy of k c \Eq. So with an increas- 
ing number of hadronic generations the particles are dis- 
tributed in more an more distinct energy levels. It turns 
out that in the generation n there are in fact n + 1 dis- 
tinct groups of particles of identical energy. This situa- 
tion is too complex to be handled in a compact analytic 
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way. According to Matthews almost all secondary par- 
ticles populate the one or two lowest energy levels per 
generation. Thus it can be justified to proceed as in the 
normal hadronic Heitler model and approximate the ter- 
minal hadronic generation n c as where the average energy 
per pion drops below the critical energy, and thus 



In ((1 + r n mult )/(l - c fv ine i)) 
where the total number of pions in generation n is 

Nl = (l + rrw t ) n 
and the average energy of these pions is 

The muon number is then 



(12) 



(13) 



(14) 



N„ 



EH 



with 
ln(l + r n mult ) 



In ((1 + r ra mu it)/ (1 - c K ine i)) 



ln(l - CK ind ) 

ln(l + rn mu i t ) 



(15) 



Since the inelasticity of hadronic interactions in air 
showers is of the order of 0.6 the approximation that the 
electromagnetic output of the first interaction is dominat- 
ing the depth of the shower maximum is in fact not justi- 
fied. Leading particles carry a large fraction of the total 
energy deeper into the atmosphere and contribute sig- 
nificantly to the electromagnetic cascade by finally gen- 
erating e.m. sub-showers at larger depths. This super- 
position of electromagnetic subshowers starting at differ- 
ent hadronic generations does not fit into the analytical 
frame of the Heitler model, and thus no relation between 
inelasticity and A" max can be given. In fact, a simple ex- 
tension of Eq. (Tin)) suggests that the value of A" max would 
decrease proportional to A ln(l — k g i), while, on the con- 
trary, it should increase with increasing elasticity because 
more energy is transferred deeper into the atmosphere. 
This is confirmed by our calculations, cf. SectionfV] In all 
these considerations we have not taken into account the 
fact that the interaction cross sections, the elasticity, the 
secondary particle multiplicity and many more parame- 
ters of hadronic interactions depend on the energy. It is 
possible to improve the Heitler model to include some of 
these effects but this goes beyond the scope of this work 
(see, for example, [6l|) . 

In Table |T] the expected relations between features of 
hadronic interactions and important air shower observ- 
ables are summarized. The results for A" max are all ob- 
tained from Eq. (|10|) . where the relation Ai oc \ja\ is 
used. The results for the muon numbers are approxima- 
tions of 



TABLE I. Proportionalities found in the extended Heitler 
model between hadronic interaction features and the aver- 
age of air shower observables. Some of the proportionalities 
are approximations that also depend on the value of other 
parameters. 
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For comparison, Matthews quotes 



that In AT,, 



oc 



« 1 - 0.14/Sinel, 

which corresponds to our result. The change of the 
electron numbers at the shower maximum is given by 

dlnAe,ma X = dN e , mox/Ne, max OC -dA^/A^ max = 



-iVA^e 



din A,, oc 



-dlnA^, which is derived from 



Eq. ©. We will confront the results listed in this table 
with our simulations in Section [V] 



D. Fluctuations of the Air Shower Development 

Not only the mean values of air shower observables are 
important, but also the fluctuations of these observables 
contain valuable information. They can complement the 
interpretation of air showers by exploiting different as- 
pects of the shower development. Compared to average 
shower properties, the investigation of shower fluctua- 
tions has the potential to yield more reliable results. 

The Heitler model has also some limited capabilities to 
explain the fluctuations of air shower cascades. Accord- 
ing to Eq. (fT0|) fluctuations in the depth of the shower 
maximum can be related to fluctuations of the depth of 
the first interaction V[Ai] = Aj oc 1/of and fluctuations 
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(16) 
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FIG. 5. Fluctuation of the longitudinal air shower develop- 
ment for 20 proton induced EAS at E = 10 19 ' 5 eV. 



6 



of the multiplicity via 

V[X max ] ocl/^ + A^Vpnn mult ]• (17) 

Thus, for small o\, the fluctuations in X max are depend- 
ing on the cross section of the cosmic ray projectile in 
the atmosphere, while for very large o\ the fluctuations 
of the multiplicity arc dominating. 

The fluctuations of the logarithm of the muon number 
can be derived from Eq. (|16p . The resulting coefficients 
arc complicated functions of c, fci ne i, n mu it and E£, and it 
is not helpful to discuss them further. Fluctuations of the 
electron number are directly linked to the fluctuations of 
the muon number and X max by Eqs. (jHJ) and Q. It is 
demonstrated in Section [V] that the correlation between 
the fluctuations in the longitudinal shower development, 
i.e. the depth of the shower maximum, and the fluctu- 
ations of the electron number is very strong and mostly 
dominates over other effects. 

To get an impression on important features of fluctua- 
tions of the longitudinal shower development of electron 
and muon number profiles, we perform some more de- 
tailed simulations beyond the Heitler model. 

The characteristic fluctuations of air showers are 
demonstrated in Fig. where electron and muon num- 
ber profiles are shown for a series of equal energy air 
showers simulated with CONEX [H. The depth of the 
electron number profiles shows large fluctuations while 
the electron number at the maximum is almost constant. 
On the other hand, the depth of the maximum in the 
muon profile docs not show large fluctuations, but the 
maximum number of muons does. Also the electron pro- 
file has a very clear maximum with a fast absorption at 
larger depth, while the muon profile only exhibits a very 
moderate absorption after the maximum. This is also re- 
flected in the logarithmic RMS of the electron and muon 
profiles around their mean profile. The electromagnetic 
part of the shower exhibit a very pronounced minimum 
in its fluctuations at the depth of the mean shower max- 
imum (X max ). 

IV. SIMULATION OF AIR SHOWERS WITH 
MODIFIED INTERACTION MODELS 

In earlier studies of similar aim it was investigated how 
the interpretation of air showers is affected by chang- 
ing internal parameters of individual hadronic interac- 
tion models |65l - [69| or by comparing the predictions of 
different interaction models j70l - t73| . These approaches 
have the virtue of not, or only marginally, leaving the al- 
lowed parameter space of the original models. However, 
underlying assumptions that cannot be justified by fun- 
damental theoretical principles, on which all models are 
built on, are often even more important than just param- 
eter values. Because of this, only the phase-space in the 
direct vicinity of the original model can be explored by 
parameter variations. It is found that typically the differ- 
ences between different models are larger than what can 



be obtained from variations of internal model parame- 
ters. Furthermore the available models do not exhaust 
the full range of current theoretical uncertainties. 

This is why for our study we are not changing inter- 
nal parameters of hadronic interaction models, but de- 
veloped an ad hoc model to explore the uncertainties be- 
yond the predictions of the available interaction models. 
Our model uses the predictions of existing hadronic event 
generators and modifies the output of these models in a 
suitable way to probe the phase-space of interaction char- 
acteristics more exhaustively. These modifications are ad 
hoc and explicitly not based on an underlying fundamen- 
tal theory or phenomenology. This allows us, firstly, to 
apply our modification to any interaction model avail- 
able and, secondly, to scan the phase-space of interaction 
physics very freely and extensively. 

It is one of the advantages of our approach that it is not 
bound to parameter correlations normally arising from 
model-internal mechanisms. Therefore we can study the 
influence of the modification of one parameter on the 
shower evolution almost independently of the other pa- 
rameters. For example, the increase of the inelasticity is 
typically correlated to an increase of the secondary parti- 
cle multiplicity. This is not the case in our way of modify- 
ing the model predictions. Of course, energy and charge 
conservation imposes some correlation between different 
parameters that are present even in our approach. 

A. Considered Parameters of Hadronic 
Interactions 

We investigate the impact of the variation of several 
important features of hadronic interactions on typical air 
shower observables. These are 

• the secondary multiplicity, rt m uit, which is defined 
as the total number of secondaries produced in a 
hadronic collision, 

• hadronic particle production cross sections, er pr od, 
that determine the interaction length of particles 
and thus the depth of the first interaction in the 
atmosphere, but also the speed of the development 
of the hadronic shower core, 

• the elasticity, K e \ = Sioading/^tot with -Evading be- 
ing the energy in the lab system carried by the lead- 
ing particle after the collision and E to t the energy 
of the projectile in the lab system, 

• and the pion charge-ratio, c = n^o / (n^o + + 
n,-), with n^x the number of pions of type 7r x , 
which determines the fraction of particles going into 
the electromagnetic cascade of the air shower after 
each interaction, assuming that all 7r° decay basi- 
cally instantly into 27. 

The air shower observables that we study are the 
shower maximum, X max , the total number of electrons 
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FIG. 6. Example of a modified cross section for Sibyll for a 
20% increase and decrease of /ig (see 23, 2j| for references 
to the data). 

above 1 MeV, N e , as well as muons above 1 GeV, N^, ar- 
riving at the slant depth of A b s = 1000 g/cm 2 , and the 
invisible energy, E- mv , which is not accessible to calori- 
mctric fluorescence telescope measurements. 

B. Extrapolation in Energy 

To investigate the importance of the extrapolation of 
high energy hadronic interaction models for the interpre- 
tation of EAS observables, we modified interaction char- 
acteristics during EAS simulations. To achieve this, in- 
dividual hadronic interaction features are altered by the 
energy-dependent factor 

f(E, f 19 ) = 1 + (/ig - 1) F{E) (18) 

with 

( E < 1 PeV 

F(E)={ log 10 (£/lPeV) E>1PcY ■ (19) 
[ log 10 (10EeV/lPeV) 

The factor f{E, f w ) is 1 below 10 15 eV, where the mod- 
els are constraint by accelerator data. Above 10 15 eV, 
the deviation of f{E, fig) from 1 increases logarithmi- 
cally with energy, reaching the value of /19 at 10 19 eV, 
cf. Fig. |6l This reflects the increasing uncertainty of 
the extrapolations with energy. The factor f(E, fig) is 
then used to re-scale specific characteristic properties of 
high energy hadronic interactions, such as the interac- 
tion cross section, secondary particle multiplicity etc.. 
Obviously, by doing this we rapidly leave the parameter 
space allowed by the original model; Thus, for large de- 
viations from the original model the results are getting 
less reliable and have to be treated with some caution. 
Nevertheless, one can get a clear impression of how the 
resulting EAS properties are depending on the specific 
hadronic interaction property. 



C. Simulation of Showers 

As primary particles we will focus on proton and iron 
at an energy of Eq = 10 195 eV. The elements in the mass 
range from proton to iron nuclei are the most abundant 
elements in cosmic rays. This way our results will bracket 
the expected mass range of the primary particles. An 
energy of about 10 19 5 eV is the highest energy at which 
one can expect to have good statistics of measured air 
showers with current and forthcoming observatories [zD, 
[75| . Above this energy the flux of cosmic rays is strongly 
suppressed pT6l. [T7j . 

The simulations are mainly based on the interaction 
model Sibyll [H, [zl, [79| . The choice of this model is 
not important for our studies as we will change the model 
output phenomenologically. It serves only as a reference 
for comparison with other models. We also performed 
cross checks with other interaction models for proton pri- 
maries. 

All shower simulations are performed with the hybrid 
code CONEX The CONEX energy threshold above 
which particles are tracked individually is set to 10 15 eV 
to require full Monte Carlo simulation for all interactions 
above 10 15 eV. The CONEX air shower simulation pro- 
gram was modified to evaluate f(E, fig) after each in- 
teraction of a hadronic particle (nucleon or meson) with 
a nucleus of air. In case f(E, fig) ^ 1 a resampling al- 
gorithm [8fjj is applied to change the properties of final 
state particles of the corresponding interaction. The re- 
sampling algorithm is designed to specifically change only 
individual properties of the secondary particle distribu- 
tions, while conserving all of the important features like 
total energy, charge, particle types, energy fractions in 
different particle types, and the leading particle as much 
as possible. For the detailed description of the resam- 
pling algorithms see Appendix [A"l 

For primary cosmic ray nuclei with mass number A > 
1, the Glauber model [HI, Ell ^ s use d to relate nucleus-air 
interactions to the underlying nucleon-air interactions. 
The hadronic event generator Sibyll is based on the 
semi-superposition model [llj], which allows us to apply 
the same resampling algorithms for showers of nuclear 
primaries. Thus, for primary cosmic ray nuclei other 
than protons we adopt the same modified Sibyll inter- 
action model but with each of the superimposed nucleon- 
air interactions being individually changed according to 
f{E, /ig)-of course at the correspondingly lower energy 
of E = Eg /A compared to the total energy of the projec- 
tile nucleus Eg. See Appendix IBI for details. 

V. RESULTS 

In the following all simulations are performed for pro- 
ton and iron primaries of 10 19 ' 5 eV. Since the results dis- 
cussed here are not very much depending on the particu- 
lar choice of the primary energy, the findings are relevant 
to the analysis of air showers at least in the energy inter- 
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FIG. 7. Impact of hadronic interaction features on the shower maximum, X max , for proton (left) and iron (right) primaries. 



val from 10 19 to 10 20 eV. For each point in the parameter 
space under investigation, 1000 showers are simulated. 

In the discussion of our results we will frequently com- 
pare to the analytic Heitler model predictions summa- 
rized in Table [I] and also refer to the dependence of EAS 
fluctuations on the longitudinal shower development as 
shown in Fig. [5J 



A. Longitudinal Shower Development and Depth 
of the Shower Maximum 

The results for the mean depth of shower maximum, 
(X max ), and the fluctuation of X max , characterized by 
RMS(X max ), are summarized in Fig. [7] The extrapola- 
tion of the total cross section for particle production has 
by far the biggest impact on X max . It can shift (X max ) 
by almost 100g/cm 2 for protons and 40g/cm 2 for iron 
in both directions, and exhibits a strong correlation with 
the fluctuations of X max . All the other interaction char- 
acteristics considered here change the fluctuations only 
within a few g/cm 2 , except the elasticity for proton pri- 



maries. A high elasticity leads to a moderate increase in 
fluctuations, at the same time shifting the (X max ) deep 
into the atmosphere. The secondary multiplicity is al- 
most as effective in shifting (A max ) as the cross section. 
This is a consequence of the distribution of the same en- 
ergy onto a growing number of particles, which is also 
predicted by the Heitler model. However, the depen- 
dence we find is somewhat different from the simple pro- 
portionality to — In n mu it for larger deviations from the 
original model. For proton primaries the dependence on 
the cross section is similar to l/cr as in the Heitler model, 
especially at larger cross sections; For iron primaries, on 
the other hand, this change is more like — In a. Further- 
more, in contrary to the independence of (X max ) from 
the pion charge ratio c we find a slight trend oc In c. The 
impact of the elasticity is approximately oc K e \. 

In addition to studying A max we also considered the 
quantity AX = A max — X\ , with X\ being the depth of 
the first interaction in a shower. AX is only sensitive to 
the shower development that follows the first interaction. 
In Fig. [5] the results for AX are summarized. 

As can be seen, only modifications of the cross section 
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FIG. 8. Impact of the modification of hadronic interaction features on AX for proton (left) and iron (right) primaries. 



change X\ and AX at the same time. By construction, 
all other modifications have an identical impact on X max 
and AX, since they do not change X\. Concerning the 
impact of the cross section extrapolation, we find that for 
iron primaries almost the full impact on X max originates 
in fact from the air shower development after the first 
interaction AX. The interaction length of iron in air 
at 10 19 5 eV is just ~ 9g/cm 2 while for protons it is ~ 
38 g/cm 2 . Reducing the cross section by a factor of /19 = 
0.5 at E = 10 19 5 eV increases the interaction length by 
a factor of 1.3, which for iron accounts only for about ~ 
12 g/cm 2 of the total effect of ~ 35 g/cm 2 seen in (X max ). 
The rest is caused by the impact of a smaller cross section 
on the air shower development, and thus AX. Even for 
proton primaries the air shower development contributes 
about half of the total effect on (X max ). For example, 
lowering the cross section by a factor of 2 increases AX 
by ~ 47 g/cm 2 while X max changes by ~ 97 g/cm 2 . 

For proton primaries the minimal possible fluctuation 
in X max are about ~ 35 g/cm 2 . This asymptotic behavior 
of RMS (Xmax) for large cross sections corresponds well to 
the predictions of the Heitler model, Eq. (|17p . where the 



asymptotic value of the fluctuations is related to fluctu- 
ations of the secondary multiplicity. For the case of iron 
primaries no such saturation effect is observed, which in- 
dicates a very much reduced importance of fluctuations 
induced by the secondary particle multiplicity. Also the 
much smaller overall dependence of RMS(X max ) on the 
cross section for iron induced showers should be noted. 

The independence of RMS(X max ) of a modification of 
the multiplicity is related to V[X max ] cx V[lnn mu i t ] (c.f. 
Eq. (JT7J)) and thus (with / being the rescaling factor) 

RMS(X max ) ocV[ln(/ ■ n mult )} = V[ln /] + V[ln n mu i t ] 
= V[lnrwt]. (20) 

A straightforward rescaling of the multiplicity, as per- 
formed by our ad hoc model, does not change the X max 
fluctuations since the relative fluctuations of the multi- 
plicity are unchanged. 

The results shown in Fig. [7] can be used to gain 
some understanding of possible interpretations of ex- 
isting data. For example, the fluctuations of X max of 
26 ± 10 (stat.) ± 5 (syst.)g/cm 2 observed by the Pierre 
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FIG. 9. Impact of a modified extrapolation of hadronic interaction features on the number of electrons, N e , for proton (left) 
and iron (right) primaries. 



Auger Collaboration [§| at 10 19 - 54 eV are compatible at 
the 1 a level with iron primaries, but also with proton pri- 
maries in combination with a strong increase of the cross 
section to higher energies. Of course, in addition to the 
I max fluctuations also the mean depth of shower maxi- 
mum and further observables have to be considered for 
obtaining a consistent interpretation of the data. Such 
an analysis is beyond the scope of this work. 

B. Electron Number at X = 1000 g/cm 2 

The impact of a modified extrapolation of hadronic in- 
teraction features on the electron number at 1000 g/cm 2 
is shown in Fig. [9] The effect on the electron number 
is twofold. Firstly, there can be some direct influence 
of the interaction characteristics on log 10 N e . Secondly, 
and generally of even more importance, the change in 
the electron number reflects also the attenuation of the 
shower in the atmosphere: when A max increases and gets 
closer to 1000 g/cm 2 the number of electrons increases, 



and when X max decreases and veers away from the obser- 
vation level, the number of electrons decreases. There- 
fore (log 10 N e ) follows very much the corresponding trend 
in (X max ), cf. Fig. [71 However, if the impact in log 10 iV e 
would be solely caused indirectly by A max then we would 
e.g. expect the fluctuations to strictly follow the trend 
shown in Fig. [5) which means fluctuations in log 10 N e are 
minimal at (X max ) and rapidly grow to smaller and larger 
depths. However, this trend is qualitatively found only 
for the modified multiplicity extrapolation. The elastic- 
ity, the charge-ratio and most clearly the cross section 
are not showing this behavior, or even an opposite effect. 
At the same time we also find that the influence of the 
multiplicity on (log 10 N c ) is largest, and the importance 
of the extrapolation of the cross section, the elasticity, 
and the charge-ratio are reduced compared to that found 
for X max . All this is evidence that in addition to the in- 
direct correlation by the longitudinal EAS development, 
the changes in the overall particle production within the 
air shower have a strong direct impact on the electron 
number. See Eq. © to get a qualitative description of 
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FIG. 10. Impact of a modified extrapolation of hadronic interaction features on the number of muons, iV^, for proton (left) 
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how -/V Cjmax and X max are both affecting the electron 
number at ground level, N e . Note, that in Fig. [9] the 
limit X bs 3> Imax is not always fulfilled. 

It is interesting to note, that an efficient way to change 
the fluctuations of the electron number in proton showers 
is via the extrapolation of the multiplicity, especially in 
the directions of reduced fluctuations. To increase the 
fluctuations, lowering the cross section can be equally 
effective. For iron primaries increasing the cross section 
leads to decreased fluctuations, while the overall impact 
of the multiplicity is very much reduced. 

C. Muon Number at X = 1000 g/cra 2 

The results of the influence of the modification of in- 
teraction features on the muon number are summarized 
in Fig. [TO] In analogy to electrons also the muon number 
at 1000 g/cm 2 reacts to changes in the depth of X max 
relative to the observation level. But as shown in Fig. [5J 
this sensitivity to X max is much smaller than in the case 



of electrons. Especially the fluctuations are not hav- 
ing the clear minimum at (X max ), but show a rather 
smooth transition to a very constant rate of fluctuations 
at larger depths. Furthermore, since muons in air show- 
ers are mainly produced via the decay of pions, their 
abundance is very sensitive to the overall number of pi- 
ons in the shower. During the shower development there 
is a competition between pion decay, yielding muons and 
neutrinos, and interaction, producing new hadronic sec- 
ondaries. While the hadronic interaction length of pi- 
ons is Aj nt ~ 120 g/cm 2 over a wide range in energy, 
the decay length changes with energy (7 = E/m) and 
is Adoc = pc7Tdcc- Muons are produced mainly by pi- 
ons with Adoc < Aint- This is, for example, why high 
energy muons are more efficiently produced at high alti- 
tude where p is small. When muons are produced, pions 
are removed from the multiplicative hadronic shower cas- 
cade and ultimately fewer pions in the air shower lead to 
the production of fewer muons. In addition, an increase 
of the production of high-energy muons at high altitude 
reduces the production of the more abundant muons at 
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lower energies. 

Since muons are hardly attenuated in the atmosphere 
(cf. Fig. [5]) the correlation of the muon number with the 
location of A max is much reduced in comparison to the 
electron number. For muons the actual particle produc- 
tion characteristics in the EAS are much more directly 
important. This is demonstrated by the large positive 
correlation of (log 10 N^) with the multiplicity, while for 
electrons a negative correlation is observed. In our sim- 
ulations the muon number increases almost proportional 
to lnrimuit. The Heitler model predicts a smaller effect 
of just — l/lnn mu it at large multiplicities. 

The muon number is also very sensitive to the charge- 
ratio. A change in the charge-ratio is directly leading 
to a shift of the fractions of the primary energy going 
into the muonic and the electromagnetic shower compo- 
nents. The Heitler model indicates log 10 oc —c and 
our results are consistent with this. 

Similar to the case for the electron numbers the pre- 
diction by the Heitler model for the dependence on the 
elasticity is log 10 oc n e \, but our simulations show 
an opposite effect more like oc — n c \. The fluctuation in 
log 10 for proton primaries depend mainly on the elas- 
ticity. 

In Fig. [11] we show the change of the muon number ob- 
tained with the modified version of Sibyll relative to the 
prediction of QGSJetII. Within the Sibyll model, the 
central parameters of the interaction have to be changed 
substantially to predict the same or even more muons 
than in QGSJetII. Only extreme modifications of the 
multiplicity or the charge ratio lead to an increase the 
muon number to the level of QGSJetII. However, one 
has to keep in mind that Sibyll is the model with the 
smallest muon numbers. 

For example Epos 1.61, which predicts the highest 
muon numbers, yields up to ~ 70 % more muons than 
Sibyll. This difference cannot be explained by a modi- 
fication of the interaction features discussed here, whose 



change would lead to an increase of the muon number 
by only ~ 30%. As noted by the authors of Epos, a 
larger number of baryon-antibaryon pairs produced in 
this model keeps a larger fraction of the primary energy 
in the multiplicative hadronic shower component to even- 
tually produce a significantly larger number of pions and 
thus muons at low energy [46j . Another important aspect 
is the chance probability of producing a neutral pion as 
a leading particle in charged pion interactions with air 
nuclei [81(. Differences in these features of hadronic in- 
teractions lead to large differences in the total number of 
muons predicted for air showers. 

Recently the Pierre Auger Collaboration reported a 
measurement of the density of muons in showers of 
10 19 eV at 1000 m from the core. Using proton induced 
showers as a reference, the number of muons is found to 
be about 1.3 times higher than that predicted by using 
QGSJetII for high energy and FLUKA for low energy 
interactions Q . Our studies show that increasing the to- 
tal muon number in showers by 30 % requires, for a given 
primary composition, rather drastic changes of the fea- 
tures of the hadronic interactions considered in this work. 
In contrast, already moderate modifications of the mod- 
elling of low energy interactions lead to a ch ang e of the 
muon number of this order, see discussion in |62l. l82l [83| . 
A good knowledge of particle production is needed both 
at high and low energy to be able to use the muon num- 
ber of showers for absolute measurements of the mass of 
the primary particle. 



D. Invisible Energy 

It is one of the advantages of the observation of air 
showers by fluorescence telescopes that the measured sig- 
nal is proportional to the energy deposit dE/dX of the 
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air shower in the atmosphere [HI, |H|, E3| • The integral 

£ c m - = JdX^(X) = JdX a cS (X)N c (X) , (21) 

where a ff is the effective (i.e. average) energy loss per 
electron, yields then an accurate measurement of the to- 
tal energy of the primary cosmic ray particle that was 
transferred to the electromagnetic shower component. 
Muons and hadrons are much less abundant compared to 
electrons and also ionize less efficiently, and neutrinos do 
not deposit any energy. The invisible energy fraction (85| 
is 

£i„v = (-Bo - E°- m -)/E and thus 

E = £ c - m 7(l - einv) (22) 

cannot be calculated in a purely model independent way, 
since e; nv has to be estimated from air shower simula- 
tions. Fortunately e; nv is only of the order of a few per- 
cent (see Fig. [TJ]) and thus the introduced model depen- 
dence of the estimated total shower energy is small. 

The invisible energy fraction is naturally correlated to 
the number of muons in the air shower. The trends we 
observe are indeed very similar to what is seen for muon 
numbers. However, the impact on e lm is generally larger 
than for N^, indicating the importance of additional ef- 
fects like enhanced muon production and differing energy 
distribution of electrons or muons. Most strikingly the 
small trend to smaller muon numbers of a rising cross 
section is reversed into a trend to larger ei nv . A large 
cross section causes an accelerated shower development. 
The rise of the invisible energy is related to an increasing 
number of high energy muons produced high up in the 
atmosphere. 

Our simulations indicate that, over the full parameter 
range of the modifications, the invisible energy fraction 
does not change by more than ~ 0.02. 



VI. DEPENDENCE ON THE INTERACTION 
MODEL 

So far all the calculations were based on Sibyll 2.1. 
We will now present a comparison of the previously 
discussed results with that obtained with the models 
QGSJetOIc, QGSJetII.3 and Epos 1.61. This study 
is limited to the case of proton primaries, since the 
Glauber model used for iron primaries does not fit into 
the framework of any of these models. 

We find that the impact of the changes of the extrapo- 
lation of hadronic multiparticlc production on air shower 
predictions is very similar for all models. Only Epos 1.61 
is found to behave somewhat different. The most striking 
difference between Epos and the other models is found 
for RMS(X max ), see Fig. [13] All model predictions seem 
to be virtually independent of the multiplicity, only EPOS 
shows a clear trend to increased shower fluctuations for 
higher multiplicities. Also the correlation of the elastic- 
ity and RMS(A" max ) is stronger in EPOS in comparison to 
the other models. Furthermore, an increase of the cross 
section within Epos does not lead to a significant reduc- 
tion of the shower fluctuations in AT max . This indicates a 
different contribution of the residual fluctuations caused 
by the secondary multiplicity in Epos than in the other 
models (c.f. Eq. p7]) ). 

The muon number itself is known to be one of the 
shower observables with the largest inherent modelling 
uncertainties. At the same time it is an excellent exam- 
ple of how model-independent the relation between mod- 
ifications of the interaction characteristics to air shower 
observables is, see Fig. [TJ] While the underlying large 
differences in the model predictions are clearly visible, 
the dependence on the modifications is very similar. 
The largest differences can be observed for the elastic- 
ity where, similar to the case of RMS(X max ), the effect 
in Epos is larger than for all other models. 
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VII. CONCLUSION AND OUTLOOK 

We developed an ad hoc model to modify features of 
hadronic interaction in the simulation of extensive air 
showers. The algorithms for these modifications were 
tested using high statistics data sets of simulated p-air 
interactions at fixed energy. We then used this ad hoc 
model to explore the impact of different extrapolations 
of hadronic interaction characteristics on the predictions 
for air shower development and typical air shower observ- 
ables. 

For simplicity it is assumed that the properties of 
hadronic interactions are sufficiently well known at en- 
ergies below 10 15 cV. Only the extrapolation of these 
features to higher energies are changed, while ensuring a 
steady transition of the description of hadronic particle 
production from low to high energy. 

Our study is mainly based on the interaction model 
Sibyll 2.1. This choice is driven by the fact that 
our model for introducing modifications to nucleus-air 
interactions is founded on the Glauber model and the 
semi-superposition model as implemented in Sibyll. 
The dependence of the results on the particular inter- 
action model is small. Cross checks with QGSJetOIc, 
QGSJetII.3, and Epos 1.61 show that predictions 
for air showers change consistently in dependence on 
the modification we introduce for hadronic interactions. 
Some deviations from a universal behaviour of the mod- 
els is found for Epos, being most likely related to the 
qualitatively different structure of this model. 

Many of the observed dependences of air shower pre- 
dictions on features of hadronic multiparticle produc- 
tion can be understood qualitatively within the simple 
Heitler model of particle cascades and its extension by 
Matthews. The most relevant results on the impact of 
different extrapolations of hadronic interaction features 
on air shower observables are 

• The longitudinal air shower fluctuations, as de- 
scribed by RMS(X max ), depend mainly on the cross 
section and less strongly on the elasticity. This 
makes fluctuations in X max a good parameter to 
study hadronic cross sections at ultra-high energies. 

• The electron number is negatively correlated with 
the multiplicity, whereas the muon number shows 
a positive correlation. The secondary particle mul- 
tiplicity provides a powerful handle to change the 
electron vs. muon number ratio in air showers. 

• The invisible energy fraction is affected by at 
most ±0.02 by uncertainties of the extrapolation 



of hadronic interaction features to air shower ener- 
gies. 

• To increase the muon number in air showers the 
multiplicity of interactions can be increased, or the 
charge-ratio decreased. A change of the predicted 
muon number in air showers of more than 30 % 
seems almost impossible to obtain by only modi- 
fying high energy interactions in the simulation as 
described here. 

It turns out that the influence of the extrapolation of 
hadronic interaction features on air shower observables 
and the dependence of these observables on the mass of 
the primary particle are of similar magnitude. Therefore 
it is very difficult to use shower measurements to esti- 
mate the characteristics of hadronic multiparticle pro- 
duction at ultra-high energy. External knowledge on the 
primary mass composition will be needed for reliable es- 
timates. Such information could come from models of 
the sources of ultra-high energy cosmic rays or the ob- 
servation of magnetic deflection of the particles in the 
Galaxy. In addition the suppression of the cosmic ray 
flux above 10 19 7 eV offers a unique window to a mono- 
or bi-clcmental cosmic ray beam: almost only proton or 
iron particles are expected to arrive at Earth [8(1 [87j at 
ultra-high energy from distant sources. 

The measurement of hadronic multiparticle production 
at LHC energy has the potential to significantly reduce 
the spread of the model divergences at cosmic ray ener- 
gies. If the information most relevant to the air shower 
development will be extracted from the LHC data, this 
will very likely lead to a breakthrough in the analysis of 
air shower data. Already now experiments at the LHC 
measured multiplicity pscudorapidity distributions in the 
central event region very accurately. However, the re- 
ally important measurements for the understanding of air 
shower cascades will be: leading particle effects, multi- 
plicities in forward direction and hadronic cross sections. 
These measurements are very difficult to perform-mostly 
because they require to measure the very forward di- 
rected particle production-and will thus take much more 
time to be measured to high precision. 
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Appendix A: Modified Characteristics of Hadronic 
Interactions 

The CONEX air shower simulation program was mod- 
ified allowing us to change specific properties of the par- 
ticle production characteristics of the underlying high en- 
ergy event generators. 

The resampling algorithms are written in order to 
specifically change individual properties of the secondary 
particle distributions, while conserving all of the impor- 
tant features like total energy, charge, particle types, en- 
ergy fractions in different particle types and the lead- 
ing particle as far as possible. Since CONEX is a ID- 
EAS simulation code, no special efforts are made to con- 
serve pj_. If the energy E of a secondary particle is al- 
tered to yi eld £ mo( j its momentum p is re-computed as 
Pmod = p V E mod -m 2 /\p\. 

The algorithms are applied to the secondary particles 
of interactions in the laboratory frame of the air shower. 
This is the frame that is most relevant for the air shower 
development, and i.e. conservation of the total energy is 
paramount within this frame. 

In the following we will describe in detail how the 
algorithms work to achieve the modification of specific 
interaction features. All algorithms are controlled by 
the energy dependent modification factor f(E, /19), see 
Eqs. (US) and JTSJ. 

For the discussion we use the quantities pseudorapid- 
ity 77 = —In tan 9/2, Ai a b = E/Eq and the number of 
charged secondaries, N^. 

1. Cross Sections 

Modifying the cross section can be implemented 
straightforwardly. No secondary particle resampling is 
needed. Only the extrapolated value of all hadronic cross 
sections has to be multiplied by f(E, fig) 

^mod = ^orgi /ig) _ (A1) 

See Fig. [5] for an example of this modification. The 
rcscaling is not only done for the primary particle, and 
hence <7 p _ a i r , but for all corresponding hadronic interac- 
tion cross sections in the EAS above the chosen transition 
energy of 10 15 eV. 

2. Secondary Multiplicity 

To change the multiplicity of secondaries, first of all, 
the leading particle needs to be excluded to prevent a 
change of the elasticity at the same time. The remaining 
particles are grouped together with respect to their type: 
nucleons, pions, kaons, photons, electrons and muons. 
To obtain a different multiplicity, particles are removed 
or duplicated at random from these groups depending 
on the aim of increasing (f(E, fig) > 1) or decreasing 



(f(E, fig) < 1) the particle multiplicity. After the par- 
ticle number has been adapted to be 

<& = r&mf 19 ) (A2) 

the kinetic energy of all particles is scaled to yield ex- 
actly the same total energy as prior to the resampling. 
Since the leading particle is spared from all these opera- 
tions, the elasticity of the interaction is conserved. Fur- 
thermore, this assures energy conservation as well as a 
constant energy ratio between the particle type groups. 
To conserve the charge, particles are changed to their 
anti-partners within particle type groups, until the total 
charge balance is restored up to a maximum charge offset 
of ±1. Finally, all particle momenta are re-computed to 
account for their modified energies. 

This minimalistic approach is optimized to conserve 
as many features of the original interaction model as 
possible. The developed method allows us to resample 
secondary particles of hadronic interactions in order to 
increase or decrease the multiplicity with as few assump- 
tions as possible. No particle energy spectra or any other 
modelling is involved. Of course artificial fluctuations can 
be introduced since only existing particles are deleted or 
duplicated. But the qualitative shape of the particle en- 
ergy spectra, the energy ratio between particle groups, 
the charge balance and the total energy as well as the 
leading particle are preserved as far as possible. 

We studied the impact of the multiplicity modification 
on high statistics simulations of proton-air collisions at 
10 19 eV. In Fig. [15] (left) it can be seen that modifica- 
tions of the multiplicity basically just scale the dA c h/d?y- 
distribution. The fact that the leading particles are pre- 
served can be clearly seen in Fig. [15] (middle) close to 
Xi a h = 1. To compensate this, there occurs an under- 
respectively over-shoot of the dA^h/dA^b-distribution 
around A] a b ~ 0.1. At smaller Aj a b the scaling of the 
multiplicity can again be identified. Modifying the mul- 
tiplicity has no impact on the resulting distribution of 

-^neutral /-^charged • 

3. Elasticity of Interactions 

To modify the elasticity k o1 = 1 - fc inol = fading /-Etot 
to obtain 

KT d = K ori S m fw) (A3) 

it is only needed to re-distribute energy between the lead- 
ing particle and the rest of the secondaries. Two facts are 
limiting the range of possible modifications of the elas- 
ticity: 

lower bound: If the total available energy of all secon- 
daries is distributed equally to all secondaries, then 
the minimal achievable elasticity, K™ od > l/n mu i t , 
of this ensemble of particles is reached. This limit 
corresponds to a complete vanishing of any leading 
particle effects. All secondaries are emitted with an 
equal share of the total energy. 
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FIG. 15. Impact of modifying the multiplicity in proton-air collision at 10 eV on the resulting distribution of secondaries. 




FIG. 16. Impact of modifying the elasticity in proton-air collision at 10 eV on the resulting distribution of secondaries. 



upper bound: Since no secondary particles should get 
deleted the energy that is bound in the mass of 
all secondary particles is not accessible to further 
increase the elasticity of the leading particle. So 
the maximal elasticity is reached when the kinetic 
energy of all particles is transferred to the leading 

Ilmult 

particle: K^ od < Ef ln /E to t- In the extreme 

i 

case this corresponds to the production of particles 
with no kinetic energy. The leading particle takes 
the full energy of the projectile minus the energy 
that went into the mass of the other secondaries. 
Since these secondaries carry no kinetic energy they 
are thus not of any relevance for the development 
of the air shower cascade. 

By design the total energy E to t is not changed during the 
procedure and since no particles are produced in addition 
to the existing ones or removed also the charge balance 
as well the particle type statistics are not altered. Be- 
cause of the changing energy of all involved particles the 
particle momenta need to get re-computed at the end. 

We studied the impact of modifying the elasticity 
on a high statistics simulation of proton-air collision 
at 10 19 cV. In Fig. [TB] (left+middlc) it is shown that a 



smaller elasticity just removes the leading particles close 
to Xi a b = 1 and slightly increases the number of parti- 
cles at smaller X^b- In the central region of the collision 
there is basically no change. On the other hand, if the 
elasticity is increased, this leads to significantly reduced 
overall particle production, since more energy is carried 
away by the leading particles. This shows also clearly 
up in the central region of the collision. The seemingly 
large effect shown in Fig. [16] (right) is due in part to us- 
ing a linear axis compared to the logarithmic axis 
in Fig. [15] (middle). 

4. Charge-Ratio of Pions 

By far most of the hadronic particle production ends 
up in secondary pions. Since there are three types of 
pions, 7T°, 7r + and ir~ , roughly a third of the produced 
particles are of each of the pion types. The charge ratio 
defined as 



with ti,x being the number of pions of type 7r x , and 
is of the order of 1/3. Since neutral pions decay al- 
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FIG. 17. Impact of modifying the charge ratio in proton-air collision at 10 eV on the resulting distribution of secondaries. 



most instantly producing 27, this is the main channel to 
transfer energy into the electromagnetic cascade of the 
shower. The charged pions, on the other side, carry on 
the hadronic shower cascade (shower core) and eventually 
decay into muons. 

We modify the charge ratio to finally yield 

c mod = c OTig /(£, / 19 ) , (A5) 

by switching neutral into charged pions (f(E, fig) < 1) 
or vice versa (f(E, fig) > 1). We apply the following 
scheme: 

positive projectiles: 7r 7r° 
negative projectiles: tt + f> ir° 
neutral projectiles: 7r ± O ir°. 

Furthermore, in all cases the leading particle is excluded 
and thus fully preserved. 

Due to this logic there are limitations in the modifica- 
tion of the charge ratio of the form 

< c mod < (A6) 
n w o + n n - + n n + 

where TT-tj-stc IS, depending on the projectile, n„+ for neg- 
ative, n n - for positive, and n n + + n n - for neutral pro- 
jectiles. Given this scheme we find the probability for 

f(E, fig) > 1: to switch one ir° — > ir STC is: 

P = {n n ° - f{E, fig) n n o)/n n o = 1 - f(E, f 19 ) 

and 

f(E, fig) < 1: to switch one 7r src -)• 7r° is: 
P= (f( E , /i9)?V - n 7r o)/ri 7r =rc . 

Since changing the pion type of individual particles 
also slightly changes their rest mass, all momenta need 
to be recomputated given the new rest mass. 

We studied the impact of modifying the charge ra- 
tio on a high statistics simulation of proton-air colli- 
sion at 10 19 eV. In Fig. Q7| (left) it is shown that the 



charge ratio just scales the dA^h/dTj-distribution, while 
in Fig. [T7] (middle) it is also clear that our strategy to 
preserve the leading particle effect works, since there is 
no impact on diV c h/dXi a b close to Xi^h = 1. Similarly 
Fig. 1171 (right) illustrates how the leading particles at high 
Xi a b are preserved and the charge ratio only changes at 
small -Xi a b- 



Appendix B: Nucleus Projectiles 

For projectiles that are nuclei of A nucleons the modifi- 
cation of interaction characteristics are performed within 
the Glauber theory [H, HU . With the semi-superposition 
model [lH it is calculated how many of the nucleons of 
the projectile nucleus, iV wounc j e d, are taking part in the 
interaction with the target nucleus (air). The interaction 
of these nucleons with the target are then computed indi- 
vidually, and can thus be rescaled by the same techniques 
as introduced in Appendix [XJ at the reduced projectile 
energy of Eq/A. 




f, 9 (p-air) 



FIG. 18. Conversion from /19 to /f 9 p using the Glauber 
theory within Sibyll. 
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The cross sections for nucleus primaries are also com- 
puted within the Glauber framework. The nucleus- 
nucleus cross section is computed based on the funda- 
mental parameters cr P ~ t p , cr P ]~ p and B c \, which are scaled 
for our purpose to 

<; p mod = m f?9 p ) <7 (bi) 
<; p m od = m av p ) <r m 

-Bela.mod = f{E, /fg P ) B c \ , (B3) 

where f p $ p is chosen to yield the following equation 

p — air/ p — p p— p o \ 

„ _ "prod I "tot, mod' °ela, mod ' -°cla,modJ 

p — air/ p — p p — p o \ ^ ' 

<rod foot , <1 > B el) 



The proton-air cross sections are all computed in the 
Glauber framework as implemented in Sibyll. This ap- 
proach assures consistency with the rest of the hadronic 
cross sections in the air shower that are modified accord- 
ing to Eq. (|A1[) . The relation between /ig and the rescal- 
ing parameter of the fundamental proton-proton interac- 
tion parameters /f 9 _p is shown in Fig. [TBI 
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